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SUMMARY 


Calculations were made of the effects of surface waviness on the external 
pressure of a supercritical airfoil at design conditions. Wave parameters 
varied include amplitude, wavelength, phase, and number of cycles. Effects of 
single and multiple waves are calculated at various chordwise locations. 
General trends of surface waviness effects on pressure distribution are 
determined and these solutions are reported. Contour deviations are imposed 
on the upper surface of the airfoil. Results are presented in a manner 
designed to facilitate ready comparison with the ideal contour static pressure 
distribution. 

INTRODUCTION 

Sharp increases in energy prices during the 1970's resulted in major 
increases in airline fuel costs. As a result, in 1976 NASA initiated the 
Aircraft Energy Efficiency (ACEE) Program. This program identified and is 
developing new technologies capable of making future transport aircraft more 
fuel efficient. 

The technology with the greatest potential for energy savings focuses on 
the reduction of viscous drag (ref. 1) by attaining laminar flow over the 
wing's surface (ref. 2). Laminar flow may occur naturally on smooth, unswept 
wings with a favorable pressure gradient. To achieve laminar flow at higher 
flight speeds, where swept wings are required, control of the boundary layer 
by suction may be necessary ( LFC ) . The feasibility of LFC was established 
during the X-21 flight program (refs. 3-8). 

Original ACEE program efforts included the design of an advanced LFC 

designated LFC 989CR9 (refs. 9-12). The design enables the attainment 
of drag divergence Mach numbers similar to that of the latest supercritical 
turbulent airfoils. At the design point, shock-free supercritical flow exists 



over about 80 percent of the upper surface and about 32 percent of the lower 
surface. The airfoil has 23° sweep, a thickness to chord ratio of 0.13, 
design Mach number of 0.755, and lift coefficient of 0.551 at a Reynolds 
number of 20 million. Features of the airfoil's design point upper surface 
static pressure distribution include a short, steep favorable gradient near 
the leading edge caused by the small leading-edge radius and a slightly 
adverse gradient over the majority of the chord with the aft region exhibiting 
a Stratford-type pressure recovery (fig. 1). The steep nose pressure gradient 
minimizes growth of cross-flow boundary layer disturbances caused by leading- 
edge sweep. The slightly adverse gradient further counteracts cross-flow 
disturbances which would otherwise require large amounts of suction to 
stabilize; also, the gradient does not severely amplify Tollmien-Schlicting 
disturbances which are more easily controlled with small amounts of suction. 
The steep aft pressure recovery is primarily designed to minimize cross-flow 
boundary layer disturbance growth (fig. 1). 

In testing this airfoil, extensive modifications to the Langley 8-Foot 
Transonic Pressure Tunnel were made. To achieve the low level of turbulence 
required for testing laminar flow airfoils, a honeycomb and five wire mesh 
screens were added to the plenum chamber. Turbulence level was reduced to 
0.05 percent. A two-wall choke was installed between the model and diffuser 
to prevent noise from propagating upstream. Finally, the tunnel walls were 
faired with a liner designed to produce an infinite two-dimensional swept wing 
flow over the model (ref. 13). 

The airfoil has a flap extending over the rearward 10 percent of the 
chord. Nearly the entire upper surface and about 85 percent of the lower 
surface is covered with individually controlled suction slots. The slots 
extend spanwise, are spaced 0.12 to 1.7 inches apart along the chord. 
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and range from 0.0025 to 0.006 inch in width. Initial objectives of the test 
include hardware verification and evaluation of the effectiveness of slot 
suction in laminarizing extensive regions of supercritical flow. 

Upon completion of the slotted model test, the upper surface of the model 
will be replaced and tested with three porous skinned panel sections fabri- 
cated by the Douglas Aircraft Company (ref. 14). The forward panel extends 
from 0 to about 20.0 inches or approximately 25 percent chord, the center 
panel extends from about 20.0 to 45.7 inches or approximately 33 percent 
chord, and the aft panel extends from about 45.7 to 69.4 inches or approxi- 
mately 31 percent chord (32% used in calculations), with the remainder of the 
chord being the flap arrangement (fig. 2). The porous skin was fabricated 
with 0.025 inch thick titanium and perforated by an electron beam to produce 
0.0025 inch diameter surface holes spaced 0.025 inch apart. A fiberglass 
structure supports the skin and blocks about 1/3 of the suction holes 
(fig. 3). The model is therefore a porous suction strip design with each 
strip approximately 1 inch wide (chord-wise). 

During fabrication of these panels, efforts were made to (1 ) build a 
surface to the closest possible tolerance, (2) obtain the smallest possible 
surface waviness, and (3) use techniques applicable to a modern aircraft 
production line. Methods for thermal stress relief, contour tolerance 
maintenance, and bonding were developed as panel fabrication progressed. 

Despite intense efforts to produce a nearly perfect porous surface, 
practical methods of building large complex LFC structural pieces resulted in 
a measurable degree of surface deviation from the ideal contour. It was 
decided therefore to conduct an analytical investigation of the effect of some 
systematic contour deviations on the external pressures in transonic flow. 

The purpose of this paper is to present these results at the airfoil’s design 
point conditions (M — 0.775, = 0.551). In the interest of providing this 
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data in a timely manner, a limited appraisal of the significance of the 
assumed contour variations is made. Upper surface results are presented as 
only that surface will be tested with porous material. 

SYMBOLS 

A wave half amplitude (see fig. 4) 

C x section lift coefficient 

Cp upper surface static pressure coefficient 

M free-stream Mach number 

X/C nondimensional chordwise location 

C chord length, “ 78.1 inches 

L half wavelength 

Acronyms 

DELTA Cp Cp with wave -Cp design 
LFC Laminar Flow Control 

PCSL percent chord starting location for wave 

RESULTS AND DISCUSSION 

This investigation of surface deviations from the theoretical LFC 989CR9 
airfoil was performed by generating sinusoidal waves and numerically super- 
imposing them onto the ideal airfoil coordinates. Wave parameters that were 
varied include half amplitude (A), half wavelength (L), phase (sine or cosine), 
and number of wave cycles (fig. 4). Table 1 gives a summary of the parameters 
used in the study, and for clarity, shows a diagram of the wave being studied. 
Data figures contain plots in ACp form; that is, the reference or design 
point ideal theoretical pressure given in figure 1 is subtracted from the 
pressure calculated at each point along the chord. This method allowed for 
greater ease in detecting small differences in the surface static pressure. 
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Calculations presented are performed by the two-dimensional analysis of 
references 15 and 16. Sinusoidal waves were generated and numerically super- 
imposed on the coordinates of the NASA LFC 989CR9 airfoil. The resulting 
geometry was input to a computer code which performs a conformal mapping of 
the airfoil using a Fourier series calculation as an initial guess after which 
a successive overrelaxation iteration technique is employed. The two- 
dimensional transonic flow is computed using potential flow equations solved 
by quadrature and finite differencing techniques of second order accuracy. 

The laminar boundary layer thickness is neglected and inviscid calculations 
yielding the static pressure distribution on a 320 point computational grid 
are obtained. The range of values of half wave amplitude and half wavelength 
in combination with wave chordwise location are selected to be representative 
of the waviness measured on the perforated wind-tunnel model panels (with 
variations). No coordinate smoothing was performed. 

Figure 5 illustrates the effect of a half cycle cosine wave of 0.005, 
0.010, and 0.020 inch half amplitude extended over an individual panel 
(forward, center, and aft). Thus, three half amplitudes (A) at three starting 
locations are given (wave location in percent of chord is given above the 
figure title). Half wavelengths vary from 20.0 to 25.7 inches, which provides 
a variation in A/L ranging from 0.00020 (0.005/25.7) to 0.00100 (0.020/20.0). 
Frame 1 (frame 1 is always to the left side of figure), with the wave over the 
panel, has A/L from 0.00025 to 0.00100; three pressure oscillations 
are produced, affecting the design C p to about 80 percent chord (where the 
sonic zone ends). Frame 2, with the wave over the center panel, has A/L 
from 0.00019 to 0.00078; the 0.00078 value produces two pressure peaks within 
the sonic zone with the largest at about 38 percent chord. Frame 3, with the 
wave over the aft panel, ranges from A/L of 0.00021 to 0.00084. The 0.020 
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inch half amplitude wave produces a single pressure peak near the aft end of 
the sonic zone; because the wave extends beyond the sonic zone, the C p 
oscillations follow. In figure 5, it is seen that for the relatively long 
wavelengths analyzed, the magnitude and nature of the resulting pressure 
variations depend not only on the surface wave amplitude but also on the 
chordwise location of the surface wave with the most severe effects occuring 
for waves near the airfoil midchord. 


Figure 6 demonstrates that the effect of a half wave of any amplitude 
that begins near the airfoil leading edge can be lessened by increasing its 
length. Frame 1 repeats the effect of a wave over the first panel. Frames 2 
and 3 indicate that the effect of a half amplitude wave as large as 0.02 inch 
becomes negligible as the half wavelength approaches the airfoil chord. 

Figure 7 shows an increasing number of cycles of a 0.02 inch amplitude 
naif cosine wave with a cycle length ranging from 25 to 33 percent chord 
(equivalent to the length of each model panel). Frame 1 repeats the 0.02 inch 
data from frame 1 of figures 5 and 6 (half wave over the forward panel) . 

Frame 2 shows that addition of a half wave over the center panel reinforces 
the pressure variations at 40 to 50 percent chord and at about 70 percent 
chord. A maximum AC p of about 0.14 occurs at 41 percent chord. Addition of 
a third cycle, equivalent to the rear panel, has no effect on the preceding 
pressure variations due to the supercritical nature of the flow but reinforces 
the previously occurring pressure variation at about 70 percent chord. These 
results illustrate the undesirable effects of successive wave cycles on 
pressure distribution. 

Figure 8 show, the effect of breaking a cosine half wave over two panels 
into two cosine half wave cycles. Two examples are presented, one for the 
forward two panels and the other for the rearward two panels. In each case. 
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the effect of the single half wave is greatly amplified by the dual half 
waves. A ACp of about 0.14 is introduced by the dual half waves at about 40 
percent chord for both cases . A second ACp peak is also introduced by the 
dual waves for each case at about 70 percent chord. From this, it can be 
deduced that shimming a model's surface to decrease wave amplitude may make 
matters worse if it has the effect of breaking the wave into more cycles 
without an appreciable reduction in wave amplitude. 

Figure 9 illustrates the effects of both sine and cosine half waves of 
0.02 inch half amplitude extended over each individual panel. Frame 1 shows 
that a sine half wave is much worse than a cosine in the nose region. Frame 2 
illustrates the sensitivity of the midchord region to both types of waves, 
although the peak AC p occurs at different chordwise locations (57 percent 
chord for the sine wave and at 38 percent chord for the cosine half wave). 

With the waves on the aft panel, the pressure variations are smaller. Figure 
9 illustrates the difficulty in predicting and understanding the complexities 
of surface pressure variations that result from surface contour deviations in 
a supercritical flow region. The maximum AC p caused by the sine or cosine 
half wave shifts, depending on wave location; panels 1 and 3 have the maximum 

ACp for the sine half wave occurring forward of the cosine half wave but 

panel 2 shows the reverse trend. 

Figure 10 is designed to show the effect of mixing half sine and half 
cosine waves and to find which combination is more detrimental to the flow at 

the design point. In frame 1, the solid line represents the 0.02 inch half 

amplitude, half cosine wave over the forward panel and the dashed line is for 
an equal amplitude sine wave over the center panel. The third line is the 
result of combining a half cosine wave over panel 1 with a half sine wave over 
panel 2 to form a full cosine wave over the total length of panels 1 and 2. 
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Note that the full cosine wave has the worst attributes of each of the 
individual half waves, but does not increase severity. Frame 2 reverses the 
location of the sine and cosine half waves. This case is more severe than 
frame 1 with the peak AC p shifted from 56 to 32 percent chord. 

Figure 1 1 compares the design point pressure distribution with and with- 
out a 0.02 inch half amplitude half cosine wave extended over the forward and 
center panels (A/L = 0.00044, PCSL = 0, see frame 2, fig. 6); and with a 2° 
downward deflection of the 10 percent chord trailing-edge flap (at constant 
lift coefficient) . Deflection of the flap decreased the effect of the surface 
wave on the chordwise pressure distribution irregularities (frame 2). Flap 
deflection, however, also changed the general chordwise pressure gradient in 
the supercritical region from adverse (positive slope) to proverse (negative 
slope). The change in pressure gradient is favorable with respect to 
Tollmein-Schlichting boundary- layer instability but probably unfavorable with 
respect to cross-flow instability. 

CONCLUDING REMARKS 

Some general trends of the effects of relatively long wavelength surface 
waviness on external pressures in the supersonic region of a supercritical 
airfoil were determined from this analytical investigation. 

(1) The magnitude and nature of the pressure variations resulting from 
surface waviness depend not only on the surface wave amplitude but also on the 
chordwise location of the surface wave with the most severe effects occurring 
for waves near the airfoil midchord. 

(2) The effect of a wave initiating near the airfoil leading edge 
decreases with an increase in wavelength. 

(3) in the airfoil nose region, an above-contour half wave has an 

appreciably greater effect on the external pressures than a below-contour half 
wave. 
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(4) The effect of multiple waves on external pressure variations is more 
severe than the effect of a single wave, with the magnitude and nature of the 
effect dependent upon the location and type of wave within the supersonic 
region. 

(5) Downward deflection of a trailing-edge flap decreased the impact of 
surface waviness and changed the general chordwise pressure gradient from 
adverse (positive slope) to proverse (negative slope). 
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A » AMPLITUDE 
L = WAVELENGTH 


LOCATION FOR WAVE 


TABLE 1 - WAVE CHARACTERISTICS SUMMARY 
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FIGURE 2 - LFC 989CR9 SUPERCRITICAL AIRFOIL, ( c~78.1 INCHES ) 



SUCTION SURFACE PANEL CONSTRUCTION 
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IGURE 3 - ELEGTRON BEAM PERFORATED TITANIUM WIND TUNNEL 
MODEL PANELS 
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FIGURE 4 - DESCRIPTION OF WAVE PARAMETERS 
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FIGURE 5- EFFECT OF COSINE WAVE LOCATION ON 
DESIGN POINT SURFACE STATIC PRESSURE. 
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FIGURE 6- EFFECT OF COSINE WAVE LENGTH AND AMPLITUDE 
ON DESIGN POINT SURFACE STATIC PRESSURE- 
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FIGURE 7. EFFECT OF NUMBER OF CYCLES OF .020 
INCH AMPLITUDE COSINE WAVE ON DESIGN POINT 
SURFACE STATIC PRESSURE- 
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COSINE/SINE SINE/COSINE 

FIGURE 10. EFFECT OF COMBINATIONS OF .020 INCH 
AMPLITUDE COSINE AND SINE WAVES ON DESIGN POINT 
SURFACE STATIC PRESSURE- 
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FIGURE 11. EFFECT OF .02 INCH AMPLITUDE , 0-58 X CHORD 
COSINE WAVE AND 2 DEGREE FLAP ON DESIGN POINT 
SURFACE STATIC PRESSURE- 




